An approximate method for treating dispersion in one-way quantum channels 
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Coupling the output of a source quantum system into a target quantum system is easily treated by 
cascaded systems theory if the intervening quantum channel is dispersionless. However, dispersion 
may be important in some transfer protocols, especially in solid-state systems. In this paper we 
show how to generalize cascaded systems theory to treat such dispersion, provided it is not too 
strong. We show that the technique also works for fermionic systems with a low flux, and can be 
extended to treat fermionic systems with large flux. To test our theory, we calculate the effect of 
dispersion on the fidelity of a simple protocol of quantum state transfer. We find good agreement 
with an approximate analytical theory that had been previously developed for this example. 
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I. INTRODUCTION 

Theoretical methods for treating non-ideal components 
in quantum networks is an important task for quantifying 
imperfections in experiments. One common example is 
photon loss in optical channels, which can be treated by 
invoking a fictitious beam splitter that mixes the chan- 
nel mode with other experimentally inaccessible modes 
. The component of the channel mode reflected by the 
beam splitter therefore corresponds to photon loss. This 
approach also allows inefficient detection to be accurately 
modeled. 

In this paper, we present a technique for treating dis- 
persion in quantum channels. Dispersion arises when 
modes acquire a phase after propagation that depends 
non-linearly on frequency. Typically, efforts are made 
to operate optical fibres at the zero-dispersion point in 
order that this effect be small, and heterogeneous struc- 
tures may be used to provide an effectively dispersion free 
channel. Nevertheless, in some circumstances, it may be 
desirable to operate in a regime where dispersion is not 
negligible. Recent proposal for implementing mesoscopic 
analogues of optical schemes, such interferometers 0,13, 
and quantum state transfer protocol Q using the quan- 
tum Hall effect will necessarily have some dispersion, due 
to the non-zero mass of quasi-clcctrons in the edge state. 
In that case, an ad hoc approach was used to estimate 
the effect of dispersion. Another quantum systems in 
which dispersion during propagation is expected to be 
important is atom lasers 

In the example of treating photon loss, an additional 
element, the beam-splitter, is added to an otherwise ideal 
channel to provide a tractable model. In analogy with 
this approach, we also introduce an extra element to an 
otherwise ideal (i.e. dispersionless) channel: a resonant, 
damped cavity operating in reflection. Near resonance, 
incident modes suffer a frequency dependent phase shift 
on reflection, depending non-linearly on their detuning 
from the resonance. This is broadly the same condition 
that arises in a dispersing channel, so the aim is to fix the 



resonance and damping of the cavity to match dispersion 
as closely as possible. 

Since there are only two parameters for the cavity, it 
is plainly not possible to treat arbitrary dispersion with 
this approach. However, we show that in simple net- 
works (without feedback or interference between differ- 
ent paths) it is possible to match up to third order in the 
dispersion relations. Thus our approach handles chan- 
nels that are not-too-dispersive, over the range of input 
frequencies. 

We begin by summarising the effects of both dispersion 
and rcfiection from a cavity. We then derive the condi- 
tions for which cavity reflection is a good approximation 
to a dispersive channel, relating the frequency and damp- 
ing of the fictitious cavity to the physical parameters de- 
scribing the dispersive channel. We then make some brief 
comments on the restrictions of this approach to channels 
in feedback systems, and fermionic systems, and derive 
a master equation for describing the dynamics for sub- 
systems connected by a one-way quantum channel. The 
paper ends with a simple example illustrating the appli- 
cation of the approach to treating quantum state transfer 
over weakly dispersive channels. 



II. PRELIMINARIES 

Consider the case of noninteracting quantum field 
propagating in one dimension. Let h = 1. Then at the 
origin (e.g. point of emission) the field can be expanded 
in terms of eigen-mode operators 

V'(0) = ^6.e— *. (1) 

UJ 

Here we are implicitly considering only modes propagat- 
ing in the positive direction. This limitation will be jus- 
tified by later (more restrictive) assumptions. The use 
of a discrete sum is for notational convenience only. A 
widely applicable expression for the dispersion relation is 
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= vk 



The group velocity is 
duj 
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where Cu is the carrier frequency. For a free nonrelativistic 
particle v ~ and a = l/2m. For an electron propagat- 
ing in an edge state typically aCii <Si so that u w t; 
(We will return later to the problem that an electron is 
not a boson.) At position L the field is 



iut+ik{uj)L 



where 



k[oj) = [2a) ^{-v + \Jv'^ + 4aw). 
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Now compare the above expressions to a dispersionless 
boson field. At the origin we again have 



</)(0)=^fe.e- 



(6) 



The (non)-dispersion relation is cj = cfc, so at position I 
the field is 



If however we also include (a) a global phase shift and (b) 
bouncing off a single-mode cavity of central frequency Wf 
and linewidth 7f then 



7f — 2i(w — Wf) 



(8) 



For this result, see for example Ref. 0. This is valid 
only if the Markovian description of the coupling of the 
external field to a single mode can be used, which requires 

Af,7f,(5w < (Ii, (9) 
where Af = (D— Wf and (5a; is the uncertainty in the energy. 

III. FEEDFORWARD 

Consider the case where the output of system s (source) 
is the input to system t (target). To model dispersion 
in the propagation between s and t we consider an non- 
dispersing reflecting off an intermediate (fictitious) cavity 
mode Cf, as shown in Fig. ^ From Eqs. Q and (jS)), this 
will work if we can make the approximation 



1 1 Z ~ _ + 6i + 2arctan-i ^ (lo) 

2a c 7f 



In this feedforward case the time delay l/c m the prop- 
agation, and the absolute phase of the field 6, are irrel- 
evant to how system t responds to the output of system 
s, as long as any classial driving fields have their timings 
and phases adjusted appropriately. Thus we can always 
choose I and 9 so that the constant and linear term in 
the expansion of the LHS of Eq. (|10|l about <D agree with 
the RHS. Thus in choosing 7f and ujf we need consider 
only higher order derivatives. Since we have two free pa- 
rameters it is natural to look at the second and third 
derivatives. Equating second and third derivative gives 

aL/u^ 



167fAf/(72+4A2)2 (11) 
167f(12Af2-7f^)/(7f^+4Af2)3 (12) 



Solving for Af and 7f yields 

= 12Af2(l + 0(V^^)), 



Af 



8La 



il + 0{^a/Lu)). 



(13) 
(14) 



The error is small when a <^ Lu. This is equivalent to 
Tp ^ Trf, where Tp = L/u is the propagation time, and 
Td = L"^ / a is the time for a pulse to disperse over a length 
scale ^ L. 

In the weak dispersion limit of 2> aui, we 
have Af/u)'^ = q{v^/aui^L) = 0{v/ujL)Oiy'^ /auj) » 
0{v/QL) = 0{l/kL). Thus from Eq. ® we have 



kL > 1. 



(15) 



(j^ In the opposite limit of <C aoj, we have A^/w^ 



0{^ya/Q/L) = 0{l/kL). Thus Eq. ^ applies in all 
regimes. It might seem surprising that our description 
puts a lower limit on the propagation distance, that it be 
much longer than a mean wavelength. This can be un- 
derstood as follows. If dispersion were significant (such 
that it is necessary to match up to the third derivative 
in Eq. (|10|l 'l over the distance of a wavelength, the prob- 
lem would be so non-Markovian that the cavity descrip- 
tion would necessarily fail. If it is deemed necessary only 
to match up to the second derivative then in principle 
Eq. H15() need not hold. However on physical grounds 
the second system cannot be within a wavelength or so 
of the first without a break-down of cascaded systems 
theory altogether. Another consideration on the limita- 
tion of validity of the theory is that for the third order 
expansion to be a good approximation we must have 

6uj < 7f. (16) 

This puts an upper bound of L which scales as {Soj)^^ . 

If all of the above conditions hold then we can write 
down a master equation for the cascaded systems s, c and 
t that will be a good description of dispersive propagation 
from s to t. 



IV. FEEDBACK OR INTERFERENCE 

In other situations the absolute time delay does matter, 
in particular with feedback. That is, if s feeds into t which 
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FIG. 1: Schematic of a triply cascaded system. The output 
of subsystem s reflects off subsystem f, and the reflected field 
drives subsystem t. No signal propagates in reverse. 



feeds back into s. In that case if wc wish to use the master 
equation description we cannot inchide a time delay l/c. 
Thus the first derivative term must come from the cavity. 
This gives 



47f 



Substituting this into Eq. Hll|) gives 



(17) 



(18) 



From Eq. Q we see that we have an inconsistency. Thus 
we cannot describe feedback for a dispersive field using 
this model. On the other hand, if a = (no disper- 
sion) then we can validly satisfy these equations with 
Af = and 7f = Av / L. Interestingly, Eq. again gives 
Eq. (O- 

Another situation where time delays matter, at least 
the difference between two time delays, is when there 
are two paths by which system s may affect system t. 
In that case, if the time difference is comparable to the 
total propagation time then the same inconsistency as 
noted above will arise. Thus the applicability of this 
approach to modelling dispersion is most promising for a 
simple forward chain, and we concentrate on this for the 
remainder of this paper. 



V. MASTER EQUATION 

To begin the quantitative analysis, we derive a general 
master equation for a triply cascaded system, shown in 
Fig. ^ where the outer systems are arbitrary, but the 
subsysteixL f plays the role of the fictitious cavity intro- 
duced to simulate dispersion. We assume that subsystem 
i £ {s, f, t} is linearly coupled to the external modes, 
according to 
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(19) 



We compute the Heisenberg equation of motion for an 
arbitrary operator, Oi of subsystem z, and make the 
Born-Markov approximation, in which we assume Ki^^ = 
■y/7i/27r is independent of uj. The resulting equation is 
a Stratonovich SDE. In order to derive a master equa- 
tion, we convert this into an Ito equation, taking care of 




FIG. 2: Fermionic dispersion treated using M-port beam 
splitters to direct modes onto separate cavities, which are 
subsequently recombined. Dotted lines represent unoccupied 
modes, and grey lines indicate weakly occupied modes. 



the spatial ordering of the three cavities (see for example 
Ref. [3). Alternatively, we can directly apply the cas- 
caded systems theory of Refs. iterating the result 
to include the third system. The master equation for the 
state matrix for the triply cascaded quantum system is 

p = -i[Hsys + H,p]+ Vi^sCs + VTfCf + VT^ctlp, (20) 
where 



H 



H.c.) 



and we have introduced the Lindblad superoperator 
I?[a]/9 = apa^ — (a^ap -I- pa^a)/2. This master equa- 
tion satisfies the requirement that dynamics in subsys- 
tem s is unaffected by the dynamics of subsystems f 
or t, and subsystem f is unaffected by subsystem t, as 
implied by the cascaded description. We have also de- 
fined the bare Hamiltonian for the uncoupled systems 
i?sys = Hs + Hf + Ht- Hs and can be arbitrary, de- 
pending on the particular application in mind. The mid- 
dle subsystem is the fictitious cavity that serves to model 
dispersion, so we take Hf = uifCfCf. 



VI. FERMIONS 

The technique described above was formulated for 
bosons. Where it breaks down for fermions is that the 
Pauli exclusion principle permits only a single particle 
per cavity mode, so that the simple linear transforma- 
tion resulting from reflection off a single cavity mode © 
does not hold. However, if there is at most one fermion 
involved in the problem, then particle statistics are irrel- 
evant and our approach can be applied. Even if there 
are many fermions, if the flux is low enough then our 
approach is applicable. Specifically, for a fermion flux 
of n per second, the average occupation of the fictitious 
cavity is at most N = n/jf, so the proposal is restricted 
to fluxes n <^ 7f. That is, n <C \/ u^/La. 

One method to extend the regime of validity of our 
method in fermionic systems is shown in Fig. [5] Here the 
output from s, plus M—l modes in the vacuum state, are 
directed through an A/-port beam splitter (Af-BS) onto 
M fletitious cavities. In this case, the average number of 
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fermions, N, is distributed over M cavities, so the mean 
occupation per cavity is N/M, which can be made small 
for sufficiently large M. The splitting is then reversed, 
and the M modes drive the final subsystem t. Physi- 
cally, it is easiest to imagine that the output of s is a 
radially symmetric mode, and that the additional AI — 1 
vacuum modes are being higher-order transverse modes. 
The fictitious A/-BS then could simply be a device that 
separates M transverse segments (e.g. wedges of a circu- 
lar wire) and sends them to M fictitious cavities. 

The procedure just described leads to the following 
master equation: 

p = -i[Hsys + H,p] 
1 ^' 

+ M 5Z + V7?Cfc + V^ct]p, (21) 

k=l 

where 



M 




k=l 



Here Hgys is as before, but with Hf = X]fc=i '^f'^l'^fe- 
might be thought that a simulation with so many sys- 
tems would be computationally expensive, but since it is 
only valid if each fictitious cavity has at most one excita- 
tion anyway, the Hilbert space dimension of the fictitious 
system as a whole is only 2*^. Moreover, the probability 
that many [that is, 0{M)] of the cavities are occupied at 
any one time is very small (since the occupation proba- 
bility N/M for any one cavity is assumed small). Thus, 
it should be possible to reduce the number of basis states 
required for a simulation dramatically. 

VII. EXAMPLE: QUANTUM STATE 
TRANSFER 

In order to demonstrate our method, we apply it to 
a proposed scheme for quantum state transfer j^] be- 
tween two remote atoms each in a separate cavity, which 
are connected by an optical channel. This scheme has 
been adapted to mesoscopic systems, using quantum dots 
instead of atoms and cavities, and quantum Hall edge 
states as a communication channel so is relevant to 
both atom-optical and solid-state systems. This system 
was sufficiently simple that it was possible to find an 
approximate analytical expression for the effect of dis- 
persion 0. Here we compare this approximation with 
the more sophisticated method we have developed here. 

The protocol works by controlling the coupling 
strength between the atom and the cavity, Qs,t{t), at 
each site in such a way that the evolution coherently 
maps excitation in one atom to excitation in the other 
atom. For an ideal channel, one class of suitable control 
pulses satisfies the relation il.s{t) = ^lt{Tp — t) = n{t). 
Dispersion in the intervening channel has two effects on 



the fidelity of the transfer protocol. Firstly, the disper- 
sion will broaden the wavepacket in the channel so that 
it will have some reduced fidelity with respect to a com- 
parable wavepacket in an ideal, dispersionless channel. 
Secondly, dispersion modifies the group velocity slightly, 
so that the wavepacket arrives at the destination at a 
slightly different time. This can be accounted for simply 
by adjusting the timing and phase of the control fields so 
that the term linear in tj — w in the expansion of Eq. H10|l . 
is zero, i.e. Tp = l/c + 47/(7^ + 4A2) k. l/c + y3/2A. 
For the purposes of feed-forward simulation, we can take 
Tp = //c = 0, so the conditions on the driving fields for 
optimal transfer is VLs{t) = rit('\/3/2A — t). 

For this model we consider Hi^s.t = ^iicjci + a|ai) -I- 
ili{t){clai + ajc;), where Ci are cavity mode annihilation 
operators, and are atomic lowering operators for each 
subsystem i, and ^i{t) is a controllable coupling between 
the atom and cavity mode. We assume the ideal case, 
= (D and "fs,t = 7- Moving to the usual interaction 
frame, the system Hamiltonian is 

Hsys = J2 + ale,) - A4c2. (22) 

i=s,t 

We assume the system starts in the state |e, 0; 0; g, 0), 
where | atoms, cavitys; cavityf; atomt, cavityt) denotes the 
states of the three subsystems expressed in the energy 
eigenbasis of the atoms and cavities. Because there is at 
most excitation, this system is equivalent to a fermion 
system 0, and there is no need for more than one ficti- 
tious cavity. 

We can now solve Eq. H2U|) for the state matrix of the 
system, which is spanned by the states 

{|g,0;0;g,0),|e,0;0;5,0),|5,l;0;5,0), 
|ff,0;l;g,0),|5,0;0;5,l),|5,0;0;e,0)}. 

We use a simple pulse sequence that implements state 
transfer n5_t(t) = 7sech(7f/2)/2 ^J. Recall that we are 
using the standard convention for cascaded systems that 
the origin of time for system t is delayed with respect to 
that for system s. 

Recall that the conditions for the cavity to accu- 
rately simulate weak dispersion are A^ = ^/Su^ /8aL and 
7^ = 12 A^, so we solve the master equation, Eq. H2()(l . 
using these parameters. In order to analyse the depen- 
dence of the infidelity, given hy = \ — T where T 
is the fidelity of the transfer, function dispersion, 
we nondimensionalise the parameters thus: a* = a/ Lu, 
A* = AL/u, 7* = jL/u. In Fig. Owe compare the re- 
sults of numerical simulations with the heuristic analytic 
expression given in 0. In that work it was found that 
the infidelity due to dispersion is given by 

J^= (a*7*^)V45, (23) 

in the weakly dispersive limit, a*7* ^ <C 1. In this regime, 
both approaches are valid and there is very good agree- 
ment, lending credibility to both. But our new method 



5 



1 -r 




6.001 0.01 0.1 1 10 100 



FIG. 3: The infidelity, 1 — versus non-dimensional diffusion 
parameter 0*7* ^. Points are from numerical calculation using 
a cavity to simulate a dispersive medium. Solid line is the 
analytic result, taken from 0. When the dispersion becomes 
dominant, the infidelity (i.e. error) asymptotes to unity. 



shows significant deviation from the approximate result 
even for a* 7*^ > 1, for which !F is still small (of or- 
der 10~^). This regime is at the limit of validity of our 
approach, according to Eq. if we say Slo ~ 7. 



VIII. CONCLUSION 

In this paper we have presented a numerical method 
for modeling the effect of dispersion in quantum chan- 
nels connecting a source system to a target system. The 
method is approximate, and can treat dispersion that is 
not too strong. We have also shown how to extend the 
approach to treat fermionic systems with large flux. Ap- 
plying our method to a simple example, for which there 
existed a previous ad hoc analytical result, showed good 
agreement between the two methods. For more compli- 
cated scenarios, analytical approaches are unlikely to be 
possible, and our technique may be the only practical 
approach. 
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